Generalizations of some integrals over the unitary 
group 



B. Schlittgen 1 and T. Wettig 1 2 

1 Department of Physics, Yale University, New Haven, CT 06520-8120, USA 

2 RIKEN-BNL Research Center, Upton, NY, 11973-5000, USA 



Abstract. Using the character expansion method, we generalize several well- 
known integrals over the unitary group to the case where general complex matrices 
appear in the integrand. These integrals are of interest in the theory of random 
matrices and may also find applications in lattice gauge theory. 
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1. Introduction 

Physical systems described by nonhcrmitian operators have recently attracted a 
lot of attention in the literature. Applications range from flux-line pinning in 
superconductors [jlj over dissipation and scattering in quantum chaos to quantum 
chromodynamics (QCD) at nonzero density In turn, the interest in these 

applications has stimulated new mathematical studies of nonhermitian matrices, in 
particular in the field of random matrix theory. 

We have been led to consider nonhermitian matrices in our work on the 
color-flavor transformation [|, ||]. This transformation involves an integration 
over complex matrices Z without any symmetry requirements. Applying the 
transformation to lattice QCD || [?], |^| results in a complex action that is not amenable 
to standard Monte-Carlo algorithms. A possible way of trying to resolve this problem 
is to write Z — HU with H hermitian and U unitary, and to perform the integral over 
U analytically so that only the integration over H has to be done numerically. This 
approach leads to the following integral over the unitary group, 

J 1= J dfi{U)det»U e^ tr< - AU+BU ^ , (1) 

U(A0 

where d/i(U) is the Haar measure of U(iV), v is an integer, which without loss of 
generality we take to be non- negative, and A,B S Gl(N, C). The above integral is the 
main focus of this paper. It is well-known for A = || [h], |ll| (where it leads to 
the effective partition function of QCD in the Leutwyler-Smilga regime [|2[ |l(|), but 
to the best of our knowledge X\ had not been computed before for A ^ B\We found 
that the latter case can be solved rather easily using the character expansion method 
put forward by Balantekin |l3|, Q . Although the case of non-integer v appears to be 
beyond this approach, we expect our result to hold in that case as well. Employing the 
same method, we could also compute the integral which arises from the factorization 
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Z = U A.V' in the color-flavor transformed action mentioned above. The matrices U 
and V are unitary, while A is diagonal with nonnegative entries. The resulting integral 
in slightly generalized form is given by 

T 2 = J dfJL(U) J dn(V) dct^UV) e^iUAVB+CVTDUi) ) (2) 
U(iV) U(JV) 

where A, B,C,D € Gl(N, C). This integral was previously known only for the case of 
v = 0, C = and D = EE, [l6|, [To). We also show that the generalization of I 2 to 
the case of unequal dimensions of U and V leads to an integral which can be non-zero 
only if determinants of U and V are absent from the integrand. We conjecture an 
expression for the result in this case. 

In Sec. ^, we derive results for T\ and X 2 . I n addition, we discuss X 2 with unequal 
dimensions of U and V and briefly consider the (trivial) generalization of the Itzykson- 
Zuber integral to the case of general complex matrices in the integrand. We also 
comment on the applicability of our results to the case of noninvertible matrices. 
Conclusions are drawn in Sec. 0. 



2. Calculation of the group integrals 

Our calculations rely on the formalism of Ref. and on several examples worked 
out in detail therein. Rather than reproducing the material from that work, we will 
refer to the relevant equations in Ref. (IJ] whenever appropriate. Thus, the reader is 
advised to have a copy of this reference at hand. 

The key observation for the generalization of the results of Ref. jL4| is that the 
representation theories of the groups U(JV) and G1(JV) are essentially the same. In 
particular, Weyl's character formula 



dct 



Xr(X) = —± L (3) 

A (Xi , . . . , XN J 

holds not just for unitary, but also for general linear matrices, cf. Ref. [|l7|. Here, 
r = (ni,...,njv) denotes an irreducible representation of Gl(iV) labelled by non- 
negative, non-increasing integers. The x t (i = X,...,N) are the eigenvalues of the 
matrix X. 



2.1. Calculation ofX\ 

Using Eq. (3.5) of Ref. 10 , we write 

det* \AU) e tr Au =YjOi r v) Xr{AU) 



tr BU' 



= 5>(°)x r (i?c/t). 



(4) 



(5) 



Here, the sums are over all irreducible representations of Gl(iV) labeled by r = 
(n%, n 2 , . . . ,ri/v) as above. The corresponding characters are denoted by Xr- F° r 
convenience, we have left out the factor of \ in the exponent of Eq. ([[]), which will be 
reinstated in the final result. The coefficients in the character expansion are given by 



A v) = det 



1 

v + i- j 



(G) 
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where i and j run from 1 to N, labelling the rows and columns of the matrix. We 
thus obtain 

2i = 



(7) 



d^U) def(AU) e t ^ AU+BU ^ 

/ dn{U)xMU)xABtf). 

f r 1 _ _ , 



U(JV) 



U(JV) 

Now, we can write 

Xr {AU)xABU^) = A%U$B^U%>, (8) 

where the superscript serves as a reminder that these matrices live in the 
representations r and r' of Gl(-ZV), respectively. If we restrict ourselves to unitary 
matrices, we have corresponding irreducible representations of the subgroup U(7V) of 
Gl(iV), and hence, we can use the group theoretical result 

1 



(9) 



U(7V) 



where d r is the dimension of r, regardless of whether we consider it as a representation 
of Gl(AT) or U(JV). It then follows that 



(o) 

Xr(AB). 



r 

From Eqs. (3.5), (3.3) and (2.9) of Ref. ju|, we obtain 

with ki = N + n.i 



(tf-i)! 



(10) 



(11) 



The matrix AB has JV eigenvalues which we denote by /Ltf , . . . , /z^y . We now apply 
Weyl's formula, 

2{ nj +N-j)~ 



det 



Xr (AB) 



Mi 



A(M 2 ) 



(12) 



where 



JY 



(13) 



is the Vandermonde determinant. Inserting this expression, together with Eqs. 
and (0), into Eq. © yields 
nv-i 



n 

n=l 



1 



E det 



%!(% - JV-i/ + i)! 



det 



Mi 



(14) 



Applying the Binet-Cauchy formula, see Eq. (B4) of Ref. [[I4[ , and using the power 
series expansion of the Bessel function, 



u 



E 

fc=0 



fc!(fe + A)! 



(15) 
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\r_i ~i 

1 



2i = 



N-l 



A(M 2 ) 



det 



(16) 



Pulling out the factors of /ij and rearranging the determinant using /„ = 7_„ yields 
det 



f jN-*+"I i _ N _ ]/ (2 lij ) =det*(AB)det 



N-i; 
L 3 



= det* (AB) det fi)' 1 (2 W ) . (17) 

Finally, we rescale A and -B by | to reinstate the factor of | in the exponent of Eq. ([!]) 
to obtain 



d M (i7) defU e |tr(AC/+BC/+) 



U(iV) 



2 — 5 — 



'JV-l 

n 



/det 5 
\deTZ 



I det 



(4 X I v +j-i{ni) 



(18) 



where, once again, the /z? are the eigenvalues of AB. 



2.2. Calculation 0/Z2 

Note first that in the integrand of I2 , the determinants of U and V have to be raised 
to the same power, as indicated in Eq. @; otherwise the integrations over the U(l) 
subgroups of U and V simply render I2 zero. 
Using again Eq. (3.5) of Ref. we have 

det" (17 AVB) e tr UAVB 



J24 u) Xr(AUVB) , 

r 

e tr OVtDfft = ^ a (0) Xr(C ryt £> f 7 t ) 
r 

with given in Eq. (^). Thus, 

i 2 = 1 y ^(l/)det-(;7 J 4yi?) e tI '( c/ ^ s + c;/tDC/t ) 

U(JV) U(2V) 

= EE«^ ) / / rfMWXr(C/AyS) Xr '(^tDC7t) 

r r ' U(IV) U(JV) 



E 



H (0) 
af a r 



d r 



U(7V) 

<z>(f/) Xr(BUADU^C) 



U(N) 



(") (0) 

= E^-^(^)^(^)> 

r r 

where we have made use of Eqs. (|) and (|). From Eq. @ we have 



(") (0) 



7V-1 

n 



-, 2 



AT 

n 



(19) 
(20) 



(21) 



(22) 
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Now denote the eigenvalues of the matrices AD and BC by x\ , . . . , x 2 N and yf, . . . , y%, 
respectively. Using again Weyl's formula (|l^), Eq. (^l|) thus becomes 

2 



x 2 = 



'N—l 

n 



n=l 



1 



En 



i 



■det 



det 



2ki 



(23) 



Employing the expansion theorem given in Eq. (B2) of Ref. [B and noting Eq. ( [l5| ) 
again, this yields 

-i 2 

L—^efct/teife)] (24) 



I 2 = 



iV-1 
n=l 



A(a: 2 )A(j/ 2 ) 



with 

f(z) = z"I- v (2z). (25) 
Pulling the factors of and y" out of the determinant and using /„ = i_ n , we obtain 

fN-l 



j 

-i 2 



n 

71=1 



det 2 (ABCD) 



det [4(2^)] 
A(x 2 )A(y 2 ) 



(26) 



We finally rescale A, £?, C and D by l/v2 to reinstate the factor of ^ in the exponent 
of Eq. (0) to obtain 

d/x(C/) / d A1 (F)det"(C/y)e^ r ^ v ' B + cv ' tl3C/t ) 



U(AT) 



U(JV) 



2iV(JV-i) 



'JV-1 
,n=l 



/ det (CD) y det [/.(a^)] 
Uet(AB) J A(z 2 )A(2, 2 ) 



(27) 



Once again, the xf and y 2 are the eigenvalues of AD and -BC, respectively. 



Let us now consider the case in which U and V have different dimensions, i.e., 



r(N,M) 



d(j,(U) 



(28) 



U(JV) U(M) 

In this case, A and C are complex N x M matrices, and B and D are complex M x N 
matrices. For definiteness, we shall take M < N. 

At this point, we have not been able to prove a result for 2^ ' > but we 
conjecture, based on explicit calculations for small N and M, as well as on numerical 
experimentation, that the result takes the form 



j-(N,M) _ 2 M ( N ~ 1 ) 



N-l 

n 



M-l 



det [I N - M (xiyj)} 



A(x 2 )A(y 2 )n^ 1 (^y,) 



N-M 



(29) 



n m! 

\_n=N-M J Lm=N-M 

Here, x\ and yf (i = 1, . . . , M) denote the (non-zero) eigenvalues of DA and BC, 
respectively. This expression also reduces to the well-known result in the case where 
C = fit and D = A\ cf. Refs. @ [HJ @. 

Note that we have not included any determinant terms in the integrand of 
2^ ■ If we included, say, det" U det n V in the integrand, integrations over the 
U(l) subgroups of U(JV) and U(M) show immediately that the value of the resulting 
integral is zero unless, possibly, v and r\ are related by Nv = Ml). We now 
show that, even if this relation holds, the integral gives zero for any v ^ 0, and 
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hence also for n 7^ 0. To see this, suppose that v 7^ 0, and let us perform the 
integral of U over U(iV), leaving the integral over U(M) untouched for the moment. 
The result could be read off from Eq. Jl8| ) if the matrices AVB and CV^D were 
Gl(iV) matrices. However, since M < N, these matrices are not of full rank, and 
therefore N — M of their eigenvalues are equal to zero. Except on a set of measure 
zero they have the same rank, so that a limiting process leads to a finite value 
of det{CV*<D)/ det(AVB), given by det(V^DC)/ det(VBA). We then need to find 

the limit of det j A(p 2 ) as (J.m+%, ■ ■ ■ , /Mr — > 0. In this context, 

jtzf , . . . , fx\n denote the non-zero eigenvalues of AVBCV^D. In fact, letting just 
fj-N — > 0, it is easy to see that the above expression goes to zero, unless v = 0, 
which shows that 2^ = 0; unless v = n = 0. 



2.3. Generalization of the Itzykson-Zuber integral 
The integral 

Is= J d/i(C/)e tr ( AC/BC/t ) = 

U(JV) 

was computed in Ref. p8[ for the case where A and B are hermitian matrices with 
real eigenvalues x\, . . . ,Xn and yi, ■ ■ ■ , yj\r, respectively. This is a special case of a 
more general result due to Harish-Chandra [jl9| . 

Following the calculation of this integral in Ref. [[u], it is immediately obvious 
that the only change in the final result is the replacement of the eigenvalues of the 
hermitian matrices A and B by the eigenvalues of their general complex versions. 

Also, including the determinant of U in the integrand (raised to any nonzero 
power) would give zero due to the integration over the U(l) subgroup. 



N-l 

n 

n=l 



det [expjxjyj)] 
A(x)A(y) 



(30) 



2.4- Comment on noninvertible matrices 

In deriving the above results, we have assumed that the matrices A, B, C and D 
are elements of Gl(N, C). However, the integrals I\ through I3 exist even if the 
matrices on which they depend are not of full rank. In this case, we can consider the 
limit in which one or more of the eigenvalues of the matrix approach zero. Uniform 
convergence permits the interchange of this limit and the integration over the unitary 
group. A PHopital procedure on the right-hand side of Eq. (18), (|27]) or ( |30| ) then 
leads to a finite (though possibly zero) result. 



3. Conclusions 



We have derived generalizations of several well-known integrals over the unitary group 
to the case where general complex matrices appear in the integrand. These integrals 
may find applications in lattice gauge theory but are also of purely mathematical 
interest, in particular in the theory of random matrices. As mentioned in the 
introduction, our motivation for studying these integrals originated from the complex 
action problem that arises if the color-flavor transformation is applied to lattice QCD. 
We found that the integral X\ solved this complex action problem for one quark flavor 
but unfortunately not for two or more flavors. 
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Our results were obtained by a straightforward application of Balantekin's 
character expansion method. It would be interesting to investigate the feasibility 
of other well-known methods to compute integrals over the unitary group, such as the 
diffusion equation method, in the cases we have considered. It should also be possible 
to generalize the results of the present paper to integrals over the super-unitary group. 
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